N93-269t>l 


Next Generation System 
Modeling of NTR Systems 

John J. Buksa and William J. Rider 
Los Alamos National Laboratory 
October 22, 1992 


- Los Alamos 


Introduction 


O NTR Modeling Challenges 

□ Current Approaches 

□ Shortcomings of Current Analysis Methods 
O Future Needs 

G Present Steps Toward These Goals 
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ENGINE COUPLING PHENOMENA 
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3D NTR Cluster 


Introduction: Modeling Applications 


• Design: performance (SS operation) and lifetime (fuel / criticality) 

• Startup and Shutdown 

(two phase T-H, neutronics, kinetics, heat transfer, low strain rate hydro) 

• Water Immersion 

(kinetics, neutronics, all hydro) 

• Impaction 

(kinetics, neutronics, high strain rate hydro) 

• Engine-Out Operations 

(all except high strain rate hydro) 
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DETAILED MCNP MODELING OF NUCLEAR THERMAL 
ROCKETS - WESTINGHOUSE NRX-A6 REACTOR 



Los Alamos 

3192 


Thermal-Hydraulic Analysis Methods 


• Extensive experience in both space and terrestrial reactors 

• TRAC 

Developed for LOCA analysis of PWRs 
Highly developed models for two-phase flow 
Low/zero gravity models are available 
Useful for facility/more general system analysis 

• HERA 

Developed for solid core terrestrial reactors 

Useful for the thermal analysis of general systems including space 
nuclear systems 

• KLAXON 

New thermal hydraulic systems code designed specifically for gas 
cooled, space reactors 

• THROHPUT 

State-of-the-art heat pipe modeling from startup to shutdown 

Los Alamos 
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• The code computes solution with minimal computatio 


Example TRAC Nodlng Diagram 








Thermal-Hydraulic Modeling: Prismatic Fuel 


• HERA: HElium/Hydrogen Reactor Analysis 

• Used to model reactor core and core components with axially homogeneous 
construction 

• Three-dimensional, fully transient, arbitrary user defined geometries 

• Programmed to be computationally efficient, especially on vector 
supercomputers 

• Currently exists in stand-alone mode and coupled to TRAC. Connection 
to KLAXON is planned 

• PATRAN grid generator and visualization translators currently being written 

• Coupling to Storm's corrosion model envisioned -► Core Lifetime 

• Component and core T-H model planned (fuel element, support element, 
and periphery) 

— Los Alamos — 


Methodology: New 
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KLAXON 

GAS-COOLED REACTOR SYSTEMS MODELING CODE 

Time-dependent analysis of systems operating with compressible gas 
working fluids. TRAC-like pipe, plenum, etc. component models, fill 
and break capabilities, and advanced flow modeling numerics for 
shock following in nozzles. 


01 

f 

Future Development | 

- Connection to HERA | oo 

- Validation with 1 

systems data | 
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THROHPUT 

HEAT PIPE MODELING CODE 


Transient thermal-hydraulic heat pipe modeling code with: 


- Multi-region capability (wall, 
fluid, mixed, gas) 


- 2-D convection and 
conduction heat transfer 

- Li melt model 

- Gravity and non-gravity 
capillary pressure models 


Future development: 
Benchmarking and validation 
with LANL experiments 


Evaporator 

Section 


Adtobattc 
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Condenser 

Section 



Heat Pipe Operation 
Los Alamos ■ 
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Future Needs 


□ Better All Around Resolution of Problems 

□ System Design Optimization Tools 

□ Complete Utilization of Modern Technology 
(Computers and Algorithms) 

□ Use of Integrated Physics Codes 
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Los Alamos Perspective 

• Emphasis on Simulation Instead of Testing 

— current ES&H environment dictates reduced testing 
of nuclear systems 


• Interagency NTP Modeling Team 

— Role, Impact, Importance, Visibility 


• Effort Should be Commensurate With the SEI 

— ambitious, high profile, high tech, national 


importance 
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Advanced Architecture: Description 


common file storage 


workstation 
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Advanced Architecture: Potential Physics Packages 


'• Neutronlcs (including cross-sections, dosimetry) 

• Spatial Kinetics 

• Generation/Depletion 

• Thermal- Hydraulics (two phase) 

• Low Strain Rate Hydro 

• High Strain Rate Hydro (solid and fluid) 

• Heat transfer (conduction, radiation) 

• Chemistry/Materlals 



• • • 
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Other Laboratory Capability 


• Fluid dynamics codes 

Developed for a large range of physical situations varying from 
incompressible to highly compressible flows 

Advanced methodologies 

• High Strain Rate Solid/Hydrodynamics 

Applicable to events involving reactor impaction/disassembly 
Examples: launch accidents, reentry, water immersion 
Coupled directly to other physical phenomena (neutronics for instance) 
Advanced methodologies 

• High Performance Computing 

One of two DOE centers of excellence 
ICN (3 CMs, 7 Cray YMPs) 

ACL 

Los Alamos — 


ADVANCED COMPUTING 
LABORATORY 

Acting as a university/ industrial/laboratory interface for 

state of the art computations, emphasizing: 

• State of the art hardware for massively parallel 
computation (largest CM-2s and CM-5 in the nation) 

• Wide area gigabit network for distributed parallel 
computing (using ANSI standard: HIPPI) 

• Advanced scientific visualization using high speed 
networking and parallel computational methods 

• Software tools/ algorithms development for distributed 
parallel computation (NSF Science & Tech, center: CRPC) 

• Emphasizing "real" applications running in parallel 
environment (Grand Challenges and beyond) 
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Purposes of the ACL 

• To respond to the rapid changes in hardware and software 

• To investigate new "Grand Challenge" computing 

environments 

• To provide more "access" to Los Alamos from the outside 

world 

• Provide high performance testbed for networking and 

visualization 

• Stimulate practical algorithm development for massively 

parallel computing 

• Function as one of the Dept of Energy High Performance 

Computing Research Centers 
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Table 1: TODAY 


HHt 

2-d immiscible 
flow 



40 GBytes 

Novel Materials 

2-d molecular 
dynamics 

wm 

500 Mbytes 

64 GBytes 


3-d multimaterial 
hydro (200 3 pts) 

wm 

8 GBytes 

lOOGBytes 

Plasma physics 

transport scaling 


8GByte8 

200 GBytes 

Global Ocean 

decade, 20 levels, 
1/2° 


500 MBytes 

250 GBytes 

Brain Topology 

3-d reconstruction 


200 MBytes 

10 GBytes 

QCD 

quenched lattice 
(32x32x32x64) 

To 15 

500 MBytes 

500 MBytes 
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Table 1: TOMORROW 


Porous Media 

3-d immiscible 
flow 


1 Tbytes 

4TByte8 

Novel Materials 

3-d molecular 
dynamics 


20 Gbytes 

3 TBytes 

3-d muKimatcrial 
hydro (1000 3 pts) 


1 TBytes 

20TBytes 

Plasma physics 

numerical Tokamak 


1 TBytes 

100 TBytes 

Global Ocean 

century, 40 levels, 
1/4° 

iHH 

4 GBytes 

20 TBytes 

Brain Topology 

3-d reconstruction 

SB 

15 GBytes 

1 TBytes 

QCD 

quenced lattice 
(64x64x64x128) 


8 GBytes 

8 TBytes 
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Applications on the CM-2 

• QCD 

• Condensed Matter Physics 

• Free Lagrange Hydrodynamics 

• Global Ocean Model 

• Lattice Gas (porous media) 

• Oil Reservoir: Mobil (llGflops sustained) 

• Tokamak Fluid Turbulence 

• Fokker Planck 

• Crystal Formation 

• Many Body Problem 

• Plasma Particle Simulations 

• Molecular Dynamics 

• Neural Networks 
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□ A 3-D Multi-Material Hydrodynamics Code on the 
Connection Machine 

□ High-Speed Hydrodynamics and High-Rate 
Deformation of Solids 

□ Euierian, Second-Order Predictor Corrector 
Lagrangian Step with Third-Order 
High-Resolution Advection 

□ High-Resolution Interface Reconstruction 
Algorithm 

□ Highly Efficient for the Connection Machine 
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The NIKE Codes 


Ptf 





D 

X2 

O 

GO 

D 

■c 8 

| g 

ie 

& S 

I s 


i 


8 ? 

Q 1 
cn u 

I I 


els 

(O eJ ’O 

I I 



Conclusions 


□ Current Modeling Approaches are Generally 
Inadequate 

G In the Future Modeling will be Relied on Heavily 

G Los Alamos has begun to Lay the Groundwork for 
Future Modeling Capabilities 
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